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ABSTRACT 


The  paper  has  two  parts.  In  part  I it  describes  .he  old  results  on  the 
problem  of  smoothing  a given  bi-infinite  sequence  of  equidistant  data  ([4],  [5], 
[1] , [2]).  It  includes  in  53  the  main  results  on  the  cardinal  spline  interpola- 
tion of  equidistant  data  of  power  growth.  All  this  is  meant  as  background  for 
the  new  developments  in  Part  II. 


In  Part  II  we  present  a method  of  smoothing  a sequence  of  equidistant  data 
{yv}  (~o°  < v < °°)  . It  is  based  on  the  ideas  of  E.  T.  Whittaker  and  may  be 
described  as  follows.  In  terms  of  the  central  B-spline  M2m(x),  defined  by  (7) 
of  Part  I,  §3,  we  define  the  cosine  polynomial  <f>2m^  bY  (26)  °£  Part  §3. 
Given  a positive  smoothing  parameter  e , we  define  the  sequence  of  "weights" 

{ ouv  ( e)  } by  the  Fourier  exoansion 
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where 


< °°  . It  is  then  shown  that  the  cardinal  spline 
^2m  * 


S(x; e)  = l 0^  (x- j ) , 


c.  = J y w.  ( e)  , 

3 L v i-v 

Vs-  -CO  J 

is  the  solution  of  the  following  minimum  problem:  Among  all  functions  f(x) 
such  that 

f(m)  (e)  e L2(3R) 

to  find  f(x)  such  that  the  functional 


(f)  = e / (x) ) dx  + y (f (v)  - y )2  is  to  be  minimal. 
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SIGNIFICANCE  AND  EXPLANATION 

Let  (y  ) (-»■  < v < “)  be  a given  bi-infinite  sequence  of  data  that  are  in 
v 

need  of  smoothing,  i.e.  replacement  by  a supposedly  smoother  sequence  y ^ . Let 
m be  a natural  number  and  e a positive  smoothing  parameter.  It  is  here  shown 
how  to  construct  a so-called  cardinal  spline  function  S(x;i)  having  the  follow- 
ing properties. 

1.  S(x*r)  reduces  to  a polynomial  of  degree  < 2m-l  in  each  of  the  inter- 
vals (i,  i+1) , for  all  integers  i . 

2.  S(x;i)  has  2m-2  continuous  derivations  for  all  real  x . 

3.  If  there  exists  a polynomial  P(x)  of  degree  < 2m-l  so  that  y^  = P(v) 

for  all  integers  v , then  S(v;e)  = y^  for  all  v . 

4.  If  e = 0 , then  S(v;0)  = y , for  all  v , provided  that  the  y do 

v v 

not  grow  too  fast  with  |v|,  in  fact  y = 0(|v|Y)  for  some  y > 0 . 

^ t\, 

5.  If  £ > 0 , then  the  smoothed  sequence  (y  1 is  given  by 

% 

y = S(v;t')  for  all  v . 


The  spline  S(x;v)  is  obtained  by  an  optimal  compromise  (depending  on  t ) 
between  two  simultaneously  irreconcilable  situations: 

(i)  That  s(x;e)  should  reduce  to  a polynomial  of  degree  £ m-1  , 

(ii)  That  S ( v; c)  = y should  hold  for  all  v . 

v 

If  we  increase  the  smoothing  parameter  t , then  we  place  more  emphasis  on  the 
condition  (i) , while  allowing  {S(v;e)}  to  differ  more  from  { y } . 

In  a future  publication  the  author  hopes  to  discuss  the  application  of 
cardinal  spline  interpolation  to  the  problem  of  approximating  the  Hilbert  trans- 
form of  a function  f(x)  defined  by  the  sequence  of  equidistant  ordinates 
{f (vh) } (-»  < v < h > 0) . If  this  sequence  is  subject  to  experimental  errors, 
it  would  seem  advisable  to  perform  on  the  sequence  { f ( vh ) ) a preliminary 
smoothing  operation  by  the  method  of  the  present  paper.  , ■ ^ 
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described  as  follows.  In  terms  of  the  central  B-spline  M^tx)  , defined  by  (7) 
of  Part  I,  §3,  we  define  the  cosine  polynomial  ^^(u)  by  (26)  Part  §3. 
Given  a positive  smoothing  parameter  r , wo  define  the  sequence  of  "weights" 

{ oj  ( c)  } by  the  Fourier  expansion 
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Let  J |y^j  < <*>  . It  is  then  shown  that  the  cardinal  spline 


S(x;  c)  = l c M^fx-j)  » 


where 


c.  = Y y a>.  ( t ) , 
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is  the  solution  of  the  following  minimum  problem:  Among  all  functions  f (x) 


such  that 

f (m)  (e)  € L2(1R) 

to  find  f (x)  such  that  the  functional 

OO 
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(f  * (x))  dx  + l (f (v)  - y ) is  to  be  minimal . 


